Structural basis and analysis of hamster ACE2 binding to different SARS-CoV-2 spike RBDs.

Pet golden hamsters were first identified being infected with the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) delta variant of concern (VOC) and transmitted the virus back to humans in Hong Kong in January 2022. Here, we studied the binding of two hamster (golden hamster and Chinese hamster) angiotensin-converting enzyme 2 (ACE2) proteins to the spike protein receptor-binding domains (RBDs) of SARS-CoV-2 prototype and eight variants, including alpha, beta, gamma, delta, and four omicron sub-variants (BA.1, BA.2, BA.3, and BA.4/BA.5). We found that the two hamster ACE2s present slightly lower affinity for the RBDs of all nine SARS-CoV-2 viruses tested than human ACE2 (hACE2). Furthermore, the similar infectivity to host cells expressing hamster ACE2s and hACE2 was confirmed with the nine pseudotyped SARS-CoV-2 viruses. Additionally, we determined two cryo-electron microscopy (EM) complex structures of golden hamster ACE2 (ghACE2)/delta RBD and ghACE2/omicron BA.3 RBD. The residues Q34 and N82, which exist in many rodent ACE2s, are responsible for the lower binding affinity of ghACE2 compared to hACE2. These findings suggest that all SARS-CoV-2 VOCs may infect hamsters, highlighting the necessity of further surveillance of SARS-CoV-2 in these animals.IMPORTANCESARS-CoV-2 can infect many domestic animals, including hamsters. There is an urgent need to understand the binding mechanism of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants to hamster receptors. Herein, we showed that two hamster angiotensin-converting enzyme 2s (ACE2s) (golden hamster ACE2 and Chinese hamster ACE2) can bind to the spike protein receptor-binding domains (RBDs) of SARS-CoV-2 prototype and eight variants and that pseudotyped SARS-CoV-2 viruses can infect hamster ACE2-expressing cells. The binding pattern of golden hamster ACE2 to SARS-CoV-2 RBDs is similar to that of Chinese hamster ACE2. The two hamster ACE2s present slightly lower affinity for the RBDs of all nine SARS-CoV-2 viruses tested than human ACE2. We solved the cryo-electron microscopy (EM) structures of golden hamster ACE2 in complex with delta RBD and omicron BA.3 RBD and found that residues Q34 and N82 are responsible for the lower binding affinity of ghACE2 compared to hACE2. Our work provides valuable information for understanding the cross-species transmission mechanism of SARS-CoV-2.

Sialoglycan binding triggers spike opening in a human coronavirus.

Coronavirus spike proteins mediate receptor binding and membrane fusion, making them prime targets for neutralizing antibodies. In the cases of severe acute respiratory syndrome coronavirus, severe acute respiratory syndrome coronavirus 2 and Middle East respiratory syndrome coronavirus, spike proteins transition freely between open and closed conformations to balance host cell attachment and immune evasion1-5. Spike opening exposes domain S1B, allowing it to bind to proteinaceous receptors6,7, and is also thought to enable protein refolding during membrane fusion4,5. However, with a single exception, the pre-fusion spike proteins of all other coronaviruses studied so far have been observed exclusively in the closed state. This raises the possibility of regulation, with spike proteins more commonly transitioning to open states in response to specific cues, rather than spontaneously. Here, using cryogenic electron microscopy and molecular dynamics simulations, we show that the spike protein of the common cold human coronavirus HKU1 undergoes local and long-range conformational changes after binding a sialoglycan-based primary receptor to domain S1A. This binding triggers the transition of S1B domains to the open state through allosteric interdomain crosstalk. Our findings provide detailed insight into coronavirus attachment, with possibilities of dual receptor usage and priming of entry as a means of immune escape.

Cryo-EM structure of SARS-CoV-2 postfusion spike in membrane.

The entry of SARS-CoV-2 into host cells depends on the refolding of the virus-encoded spike protein from a prefusion conformation, which is metastable after cleavage, to a lower-energy stable postfusion conformation1,2. This transition overcomes kinetic barriers for fusion of viral and target cell membranes3,4. Here we report a cryogenic electron microscopy (cryo-EM) structure of the intact postfusion spike in a lipid bilayer that represents the single-membrane product of the fusion reaction. The structure provides structural definition of the functionally critical membrane-interacting segments, including the fusion peptide and transmembrane anchor. The internal fusion peptide forms a hairpin-like wedge that spans almost the entire lipid bilayer and the transmembrane segment wraps around the fusion peptide at the last stage of membrane fusion. These results advance our understanding of the spike protein in a membrane environment and may guide development of intervention strategies.

Cryo-EM reveals binding of linoleic acid to SARS-CoV-2 spike glycoprotein, suggesting an antiviral treatment strategy.

The COVID-19 pandemic and concomitant lockdowns presented a global health challenge and triggered unprecedented research efforts to elucidate the molecular mechanisms and pathogenicity of SARS-CoV-2. The spike glycoprotein decorating the surface of SARS-CoV-2 virions is a prime target for vaccine development, antibody therapy and serology as it binds the host cell receptor and is central for viral cell entry. The electron cryo-microscopy structure of the spike protein revealed a hydrophobic pocket in the receptor-binding domain that is occupied by an essential fatty acid, linoleic acid (LA). The LA-bound spike protein adopts a non-infectious locked conformation which is more stable than the infectious form and shields important immunogenic epitopes. Here, the impact of LA binding on viral infectivity and replication, and the evolutionary conservation of the pocket in other highly pathogenic coronaviruses, including SARS-CoV-2 variants of concern (VOCs), are reviewed. The importance of LA metabolic products, the eicosanoids, in regulating the human immune response and inflammation is highlighted. Lipid and fatty-acid binding to a hydrophobic pocket in proteins on the virion surface appears to be a broader strategy employed by viruses, including picornaviruses and Zika virus. Ligand binding stabilizes their protein structure and assembly, and downregulates infectivity. In the case of rhinoviruses, this has been exploited to develop small-molecule antiviral drugs that bind to the hydrophobic pocket. The results suggest a COVID-19 antiviral treatment based on the LA-binding pocket.

Cryo-EM structure of a SARS-CoV-2 omicron spike protein ectodomain.

The omicron variant of SARS-CoV-2 has been spreading rapidly across the globe. The virus-surface spike protein plays a critical role in the cell entry and immune evasion of SARS-CoV-2. Here we determined the 3.0 Å cryo-EM structure of the omicron spike protein ectodomain. In contrast to the original strain of SARS-CoV-2 where the receptor-binding domain (RBD) of the spike protein takes a mixture of open ("standing up") and closed ("lying down") conformations, the omicron spike molecules are predominantly in the open conformation, with one upright RBD ready for receptor binding. The open conformation of the omicron spike is stabilized by enhanced inter-domain and inter-subunit packing, which involves new mutations in the omicron strain. Moreover, the omicron spike has undergone extensive mutations in RBD regions where known neutralizing antibodies target, allowing the omicron variant to escape immune surveillance aimed at the original viral strain. The stable open conformation of the omicron spike sheds light on the cell entry and immune evasion mechanisms of the omicron variant.

Conformational Flexibility and Local Frustration in the Functional States of the SARS-CoV-2 Spike B.1.1.7 and B.1.351 Variants: Mutation-Induced Allosteric Modulation Mechanism of Functional Dynamics and Protein Stability.

Structural and functional studies of the SARS-CoV-2 spike proteins have recently determined distinct functional states of the B.1.1.7 and B.1.351 spike variants, providing a molecular framework for understanding the mechanisms that link the effect of mutations with the enhanced virus infectivity and transmissibility. A detailed dynamic and energetic analysis of these variants was undertaken in the present work to quantify the effects of different mutations on functional conformational changes and stability of the SARS-CoV-2 spike protein. We employed the efficient and accurate coarse-grained (CG) simulations of multiple functional states of the D614G mutant, B.1.1.7 and B.1.351 spike variants to characterize conformational dynamics of the SARS-CoV-2 spike proteins and identify dynamic signatures of the functional regions that regulate transitions between the closed and open forms. By combining molecular simulations with full atomistic reconstruction of the trajectories and the ensemble-based mutational frustration analysis, we characterized how the intrinsic flexibility of specific spike regions can control functional conformational changes required for binding with the host-cell receptor. Using the residue-based mutational scanning of protein stability, we determined protein stability hotspots and identified potential energetic drivers favoring the receptor-accessible open spike states for the B.1.1.7 and B.1.351 spike variants. The results suggested that modulation of the energetic frustration at the inter-protomer interfaces can serve as a mechanism for allosteric couplings between mutational sites and the inter-protomer hinges of functional motions. The proposed mechanism of mutation-induced energetic frustration may result in greater adaptability and the emergence of multiple conformational states in the open form. This study suggested that SARS-CoV-2 B.1.1.7 and B.1.351 variants may leverage the intrinsic plasticity of functional regions in the spike protein for mutation-induced modulation of protein dynamics and allosteric regulation to control binding with the host cell receptor.

Site Density Functional Theory and Structural Bioinformatics Analysis of the SARS-CoV Spike Protein and hACE2 Complex.

The entry of the SARS-CoV-2, a causative agent of COVID-19, into human host cells is mediated by the SARS-CoV-2 spike (S) glycoprotein, which critically depends on the formation of complexes involving the spike protein receptor-binding domain (RBD) and the human cellular membrane receptor angiotensin-converting enzyme 2 (hACE2). Using classical site density functional theory (SDFT) and structural bioinformatics methods, we investigate binding and conformational properties of these complexes and study the overlooked role of water-mediated interactions. Analysis of the three-dimensional reference interaction site model (3DRISM) of SDFT indicates that water mediated interactions in the form of additional water bridges strongly increases the binding between SARS-CoV-2 spike protein and hACE2 compared to SARS-CoV-1-hACE2 complex. By analyzing structures of SARS-CoV-2 and SARS-CoV-1, we find that the homotrimer SARS-CoV-2 S receptor-binding domain (RBD) has expanded in size, indicating large conformational change relative to SARS-CoV-1 S protein. Protomer with the up-conformational form of RBD, which binds with hACE2, exhibits stronger intermolecular interactions at the RBD-ACE2 interface, with differential distributions and the inclusion of specific H-bonds in the CoV-2 complex. Further interface analysis has shown that interfacial water promotes and stabilizes the formation of CoV-2/hACE2 complex. This interaction causes a significant structural rigidification of the spike protein, favoring proteolytic processing of the S protein for the fusion of the viral and cellular membrane. Moreover, conformational dynamics simulations of RBD motions in SARS-CoV-2 and SARS-CoV-1 point to the role in modification of the RBD dynamics and their impact on infectivity.

Receptor binding and complex structures of human ACE2 to spike RBD from omicron and delta SARS-CoV-2.

The coronavirus disease 2019 (COVID-19) pandemic continues worldwide with many variants arising, some of which are variants of concern (VOCs). A recent VOC, omicron (B.1.1.529), which obtains a large number of mutations in the receptor-binding domain (RBD) of the spike protein, has risen to intense scientific and public attention. Here, we studied the binding properties between the human receptor ACE2 (hACE2) and the VOC RBDs and resolved the crystal and cryoelectron microscopy structures of the omicron RBD-hACE2 complex as well as the crystal structure of the delta RBD-hACE2 complex. We found that, unlike alpha, beta, and gamma, omicron RBD binds to hACE2 at a similar affinity to that of the prototype RBD, which might be due to compensation of multiple mutations for both immune escape and transmissibility. The complex structures of omicron RBD-hACE2 and delta RBD-hACE2 reveal the structural basis of how RBD-specific mutations bind to hACE2.

E484K and N501Y SARS-CoV 2 spike mutants Increase ACE2 recognition but reduce affinity for neutralizing antibody.

SARS-CoV2 mutants B.1.1.7, B.1.351, and P.1 contain a key mutation N501Y. B.1.135 and P.1 lineages have another mutation, E484K. Here, we decode the effect of these two mutations on the host receptor, ACE2, and neutralizing antibody (B38) recognition. The N501Y RBD mutant binds to ACE2 with higher affinity due to improved π-π stacking and π-cation interactions. The higher binding affinity of the E484K mutant is caused due to the formation of additional hydrogen bond and salt-bridge interactions with ACE2. Both the mutants bind to the B38 antibody with reduced affinity due to the loss of several hydrogen-bonding interactions. The insights obtained from the study are crucial to interpret the increased transmissibility and reduced neutralization efficacy of rapidly emerging SARS-CoV2 VOCs.

Cryo-EM structure determination of small proteins by nanobody-binding scaffolds (Legobodies).

We describe a general method that allows structure determination of small proteins by single-particle cryo-electron microscopy (cryo-EM). The method is based on the availability of a target-binding nanobody, which is then rigidly attached to two scaffolds: 1) a Fab fragment of an antibody directed against the nanobody and 2) a nanobody-binding protein A fragment fused to maltose binding protein and Fab-binding domains. The overall ensemble of ∼120 kDa, called Legobody, does not perturb the nanobody-target interaction, is easily recognizable in EM images due to its unique shape, and facilitates particle alignment in cryo-EM image processing. The utility of the method is demonstrated for the KDEL receptor, a 23-kDa membrane protein, resulting in a map at 3.2-Šoverall resolution with density sufficient for de novo model building, and for the 22-kDa receptor-binding domain (RBD) of SARS-CoV-2 spike protein, resulting in a map at 3.6-Šresolution that allows analysis of the binding interface to the nanobody. The Legobody approach thus overcomes the current size limitations of cryo-EM analysis.

A highly potent antibody effective against SARS-CoV-2 variants of concern.

Control of the ongoing SARS-CoV-2 pandemic is endangered by the emergence of viral variants with increased transmission efficiency, resistance to marketed therapeutic antibodies, and reduced sensitivity to vaccine-induced immunity. Here, we screen B cells from COVID-19 donors and identify P5C3, a highly potent and broadly neutralizing monoclonal antibody with picomolar neutralizing activity against all SARS-CoV-2 variants of concern (VOCs) identified to date. Structural characterization of P5C3 Fab in complex with the spike demonstrates a neutralizing activity defined by a large buried surface area, highly overlapping with the receptor-binding domain (RBD) surface necessary for ACE2 interaction. We further demonstrate that P5C3 shows complete prophylactic protection in the SARS-CoV-2-infected hamster challenge model. These results indicate that P5C3 opens exciting perspectives either as a prophylactic agent in immunocompromised individuals with poor response to vaccination or as combination therapy in SARS-CoV-2-infected individuals.

Molecular insights into receptor binding of recent emerging SARS-CoV-2 variants.

Multiple SARS-CoV-2 variants of concern (VOCs) have been emerging and some have been linked to an increase in case numbers globally. However, there is yet a lack of understanding of the molecular basis for the interactions between the human ACE2 (hACE2) receptor and these VOCs. Here we examined several VOCs including Alpha, Beta, and Gamma, and demonstrate that five variants receptor-binding domain (RBD) increased binding affinity for hACE2, and four variants pseudoviruses increased entry into susceptible cells. Crystal structures of hACE2-RBD complexes help identify the key residues facilitating changes in hACE2 binding affinity. Additionally, soluble hACE2 protein efficiently prevent most of the variants pseudoviruses. Our findings provide important molecular information and may help the development of novel therapeutic and prophylactic agents targeting these emerging mutants.

Potential Molecular Mechanisms of Rare Anti-Tumor Immune Response by SARS-CoV-2 in Isolated Cases of Lymphomas.

Recently, two cases of complete remission of classical Hodgkin lymphoma (cHL) and follicular lymphoma (FL) after SARS-CoV-2 infection were reported. However, the precise molecular mechanism of this rare event is yet to be understood. Here, we hypothesize a potential anti-tumor immune response of SARS-CoV-2 and based on a computational approach show that: (i) SARS-CoV-2 Spike-RBD may bind to the extracellular domains of CD15, CD27, CD45, and CD152 receptors of cHL or FL and may directly inhibit cell proliferation. (ii) Alternately, upon internalization after binding to these CD molecules, the SARS-CoV-2 membrane (M) protein and ORF3a may bind to gamma-tubulin complex component 3 (GCP3) at its tubulin gamma-1 chain (TUBG1) binding site. (iii) The M protein may also interact with TUBG1, blocking its binding to GCP3. (iv) Both the M and ORF3a proteins may render the GCP2-GCP3 lateral binding where the M protein possibly interacts with GCP2 at its GCP3 binding site and the ORF3a protein to GCP3 at its GCP2 interacting residues. (v) Interactions of the M and ORF3a proteins with these gamma-tubulin ring complex components potentially block the initial process of microtubule nucleation, leading to cell-cycle arrest and apoptosis. (vi) The Spike-RBD may also interact with and block PD-1 signaling similar to pembrolizumab and nivolumab- like monoclonal antibodies and may induce B-cell apoptosis and remission. (vii) Finally, the TRADD interacting "PVQLSY" motif of Epstein-Barr virus LMP-1, that is responsible for NF-kB mediated oncogenesis, potentially interacts with SARS-CoV-2 Mpro, NSP7, NSP10, and spike (S) proteins, and may inhibit the LMP-1 mediated cell proliferation. Taken together, our results suggest a possible therapeutic potential of SARS-CoV-2 in lymphoproliferative disorders.

Possible Link between Higher Transmissibility of Alpha, Kappa and Delta Variants of SARS-CoV-2 and Increased Structural Stability of Its Spike Protein and hACE2 Affinity.

The Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) outbreak in December 2019 has caused a global pandemic. The rapid mutation rate in the virus has created alarming situations worldwide and is being attributed to the false negativity in RT-PCR tests. It has also increased the chances of reinfection and immune escape. Recently various lineages namely, B.1.1.7 (Alpha), B.1.617.1 (Kappa), B.1.617.2 (Delta) and B.1.617.3 have caused rapid infection around the globe. To understand the biophysical perspective, we have performed molecular dynamic simulations of four different spikes (receptor binding domain)-hACE2 complexes, namely wildtype (WT), Alpha variant (N501Y spike mutant), Kappa (L452R, E484Q) and Delta (L452R, T478K), and compared their dynamics, binding energy and molecular interactions. Our results show that mutation has caused significant increase in the binding energy between the spike and hACE2 in Alpha and Kappa variants. In the case of Kappa and Delta variants, the mutations at L452R, T478K and E484Q increased the stability and intra-chain interactions in the spike protein, which may change the interaction ability of neutralizing antibodies to these spike variants. Further, we found that the Alpha variant had increased hydrogen interaction with Lys353 of hACE2 and more binding affinity in comparison to WT. The current study provides the biophysical basis for understanding the molecular mechanism and rationale behind the increase in the transmissivity and infectivity of the mutants compared to wild-type SARS-CoV-2.

Structural Insight into the Binding of Cyanovirin-N with the Spike Glycoprotein, Mpro and PLpro of SARS-CoV-2: Protein-Protein Interactions, Dynamics Simulations and Free Energy Calculations.

The emergence of COVID-19 continues to pose severe threats to global public health. The pandemic has infected over 171 million people and claimed more than 3.5 million lives to date. We investigated the binding potential of antiviral cyanobacterial proteins including cyanovirin-N, scytovirin and phycocyanin with fundamental proteins involved in attachment and replication of SARS-CoV-2. Cyanovirin-N displayed the highest binding energy scores (-16.8 ± 0.02 kcal/mol, -12.3 ± 0.03 kcal/mol and -13.4 ± 0.02 kcal/mol, respectively) with the spike protein, the main protease (Mpro) and the papainlike protease (PLpro) of SARS-CoV-2. Cyanovirin-N was observed to interact with the crucial residues involved in the attachment of the human ACE2 receptor. Analysis of the binding affinities calculated employing the molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) approach revealed that all forms of energy, except the polar solvation energy, favourably contributed to the interactions of cyanovirin-N with the viral proteins. With particular emphasis on cyanovirin-N, the current work presents evidence for the potential inhibition of SARS-CoV-2 by cyanobacterial proteins, and offers the opportunity for in vitro and in vivo experiments to deploy the cyanobacterial proteins as valuable therapeutics against COVID-19.

Cryo-EM density maps adjustment for subtraction, consensus and sharpening.

Electron cryomicroscopy (cryo-EM) has emerged as a powerful structural biology instrument to solve near-atomic three-dimensional structures. Despite the fast growth in the number of density maps generated from cryo-EM data, comparison tools among these reconstructions are still lacking. Current proposals to compare cryo-EM data derived volumes perform map subtraction based on adjustment of each volume grey level to the same scale. We present here a more sophisticated way of adjusting the volumes before comparing, which implies adjustment of grey level scale and spectrum energy, but keeping phases intact inside a mask and imposing the results to be strictly positive. The adjustment that we propose leaves the volumes in the same numeric frame, allowing to perform operations among the adjusted volumes in a more reliable way. This adjustment can be a preliminary step for several applications such as comparison through subtraction, map sharpening, or combination of volumes through a consensus that selects the best resolved parts of each input map. Our development might also be used as a sharpening method using an atomic model as a reference. We illustrate the applicability of this algorithm with the reconstructions derived of several experimental examples. This algorithm is implemented in Xmipp software package and its applications are user-friendly accessible through the cryo-EM image processing framework Scipion.

Site-Specific Steric Control of SARS-CoV-2 Spike Glycosylation.

A central tenet in the design of vaccines is the display of native-like antigens in the elicitation of protective immunity. The abundance of N-linked glycans across the SARS-CoV-2 spike protein is a potential source of heterogeneity among the many different vaccine candidates under investigation. Here, we investigate the glycosylation of recombinant SARS-CoV-2 spike proteins from five different laboratories and compare them against S protein from infectious virus, cultured in Vero cells. We find patterns that are conserved across all samples, and this can be associated with site-specific stalling of glycan maturation that acts as a highly sensitive reporter of protein structure. Molecular dynamics simulations of a fully glycosylated spike support a model of steric restrictions that shape enzymatic processing of the glycans. These results suggest that recombinant spike-based SARS-CoV-2 immunogen glycosylation reproducibly recapitulates signatures of viral glycosylation.

Bioinformatics Analysis of Spike Proteins of Porcine Enteric Coronaviruses.

This article is aimed at analyzing the structure and function of the spike (S) proteins of porcine enteric coronaviruses, including transmissible gastroenteritis virus (TGEV), porcine epidemic diarrhea virus (PEDV), porcine deltacoronavirus (PDCoV), and swine acute diarrhea syndrome coronavirus (SADS-CoV) by applying bioinformatics methods. The physical and chemical properties, hydrophilicity and hydrophobicity, transmembrane region, signal peptide, phosphorylation and glycosylation sites, epitope, functional domains, and motifs of S proteins of porcine enteric coronaviruses were predicted and analyzed through online software. The results showed that S proteins of TGEV, PEDV, SADS-CoV, and PDCoV all contained transmembrane regions and signal peptide. TGEV S protein contained 139 phosphorylation sites, 24 glycosylation sites, and 53 epitopes. PEDV S protein had 143 phosphorylation sites, 22 glycosylation sites, and 51 epitopes. SADS-CoV S protein had 109 phosphorylation sites, 20 glycosylation sites, and 43 epitopes. PDCoV S protein had 124 phosphorylation sites, 18 glycosylation sites, and 52 epitopes. Moreover, TGEV, PEDV, and PDCoV S proteins all contained two functional domains and two motifs, spike_rec_binding and corona_S2. The corona_S2 consisted of S2 subunit heptad repeat 1 (HR1) and S2 subunit heptad repeat 2 (HR2) region profiles. Additionally, SADS-CoV S protein was predicted to contain only one functional domain, the corona_S2. This analysis of the biological functions of porcine enteric coronavirus spike proteins can provide a theoretical basis for the design of antiviral drugs.

Mutation in a SARS-CoV-2 Haplotype from Sub-Antarctic Chile Reveals New Insights into the Spike's Dynamics.

The emergence of SARS-CoV-2 variants, as observed with the D614G spike protein mutant and, more recently, with B.1.1.7 (501Y.V1), B.1.351 (501Y.V2) and B.1.1.28.1 (P.1) lineages, represent a continuous threat and might lead to strains of higher infectivity and/or virulence. We report on the occurrence of a SARS-CoV-2 haplotype with nine mutations including D614G/T307I double-mutation of the spike. This variant expanded and completely replaced previous lineages within a short period in the subantarctic Magallanes Region, southern Chile. The rapid lineage shift was accompanied by a significant increase of cases, resulting in one of the highest incidence rates worldwide. Comparative coarse-grained molecular dynamic simulations indicated that T307I and D614G belong to a previously unrecognized dynamic domain, interfering with the mobility of the receptor binding domain of the spike. The T307I mutation showed a synergistic effect with the D614G. Continuous surveillance of new mutations and molecular analyses of such variations are important tools to understand the molecular mechanisms defining infectivity and virulence of current and future SARS-CoV-2 strains.

Understanding structural malleability of the SARS-CoV-2 proteins and relation to the comorbidities.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), a causative agent of the coronavirus disease (COVID-19), is a part of the $\beta $-Coronaviridae family. The virus contains five major protein classes viz., four structural proteins [nucleocapsid (N), membrane (M), envelop (E) and spike glycoprotein (S)] and replicase polyproteins (R), synthesized as two polyproteins (ORF1a and ORF1ab). Due to the severity of the pandemic, most of the SARS-CoV-2-related research are focused on finding therapeutic solutions. However, studies on the sequences and structure space throughout the evolutionary time frame of viral proteins are limited. Besides, the structural malleability of viral proteins can be directly or indirectly associated with the dysfunctionality of the host cell proteins. This dysfunctionality may lead to comorbidities during the infection and may continue at the post-infection stage. In this regard, we conduct the evolutionary sequence-structure analysis of the viral proteins to evaluate their malleability. Subsequently, intrinsic disorder propensities of these viral proteins have been studied to confirm that the short intrinsically disordered regions play an important role in enhancing the likelihood of the host proteins interacting with the viral proteins. These interactions may result in molecular dysfunctionality, finally leading to different diseases. Based on the host cell proteins, the diseases are divided in two distinct classes: (i) proteins, directly associated with the set of diseases while showing similar activities, and (ii) cytokine storm-mediated pro-inflammation (e.g. acute respiratory distress syndrome, malignancies) and neuroinflammation (e.g. neurodegenerative and neuropsychiatric diseases). Finally, the study unveils that males and postmenopausal females can be more vulnerable to SARS-CoV-2 infection due to the androgen-mediated protein transmembrane serine protease 2.